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The  opposed  jet  configuration  presents  an  attractive  canonical  geometry  for  the  evaluation  of  burning 
properties  of  turbulent  flames  with  past  studies  typically  limited  to  low  Reynolds  numbers.  Fractal  grid 
generated  turbulence  was  used  to  remove  the  low  turbulence  level  limitations  associated  with  conven¬ 
tional  perforated  plate  generators  with  the  turbulent  Reynolds  number  range  moved  from  50-120  to 
130-318.  Optimal  grid  configurations  were  determined  with  particular  emphasis  on  reducing  the  impact 
of  the  flow  upstream  of  the  turbulence  generators  in  order  to  facilitate  simpler  boundary  conditions  for 
computational  studies.  The  resulting  flow  structures  were  analysed  using  proper  orthogonal  decomposi¬ 
tion  and  conditional  proper  orthogonal  decomposition.  Velocity  and  reaction  progress  variable  statistics, 
including  conditional  velocities  and  scalar  fluxes,  are  reported  for  fuel  lean  methane,  ethylene  and  pro¬ 
pane  flames  approaching  extinction.  The  instrumentation  comprised  particle  image  velocimetry  with 
the  flows  to  both  nozzles  seeded  with  1  pm  silicon  oil  droplets  or  3  pm  A1203  particles.  Probability  den¬ 
sity  functions  were  determined  for  the  instantaneous  location  of  the  stagnation  point  to  eliminate  the 
possibility  of  low  frequency  bulk  motion  distorting  velocity  statistics.  Probability  density  functions  of 
flame  curvature  were  determined  using  a  multi-step  flame  front  detection  algorithm  with  estimates  of 
the  turbulent  burning  velocity  provided  along  with  a  discussion  of  alternative  determination  methods. 
The  data  sets  show  that  fractal  grids  generate  multi-scale  broadband  turbulence  and  present  an  oppor¬ 
tunity  for  a  systematic  evaluation  of  calculation  methods  for  premixed  turbulent  flames  that  undergo 
a  transition  from  non-gradient  to  gradient  turbulent  transport  while  approaching  extinction. 
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1.  Introduction 

Opposed  jet  geometries  have  been  used  extensively  to  investi¬ 
gate  premixed  and  non-premixed  combustion  under  laminar  and 
turbulent  flow  conditions.  The  comparatively  simple  flame  stabili¬ 
sation  method  and  essentially  adiabatic  conditions  lead  to  flame 
dynamics  and  extinction  being  related  to  the  aerothermochemistry 
of  the  combustion  process  rather  than  heat  losses.  Combined  with 
the  comparatively  simple  boundary  conditions,  the  opposed  jet 
geometry  provides  an  attractive  standard  test  case  for  the 
assessment  of  fuel  effects  and  closure  approximations  as  proposed 
by  Bray  et  al.  [1],  Lindstedt  and  Vaos  [2,3]  and  subsequently  by 
Geyer  et  al.  [4,5]  in  the  context  of  large  eddy  simulations  (LES). 
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Preparatory  isothermal  flow  studies  have  also  been  presented  by 
a  number  of  investigators  (e.g.,  Geipel  et  al.  [6],  Kostiuk  et  al.  [7], 
Korusoy  and  Whitelaw  [8],  and  Lindstedt  et  al.  [9])  and  flame 
structure  studies  include  velocity  and  scalar  field  measurements 
by  MounaYm-Rouselle  and  Gokalp  [10],  Kostiuk  et  al.  [1 1  ]  and  Lind¬ 
stedt  et  al.  [12],  Strain  effects  and  extinction  and  relight  character¬ 
istics  were  discussed  by  Kostiuk  et  al.  [13],  Mounaim-Rouselle  and 
Gokalp  [14],  Sardi  and  Whitelaw  15],  Korusoy  and  Whitelaw  [16] 
and  Luff  et  al.  [17]  amongst  others.  Mounaim-Rouselle  and  Gokalp 
[10],  Kostiuk  et  al.  [13]  and  Kostiuk  [18]  showed  that  larger  nozzle 
separations  caused  a  low  frequency  axial  movement  made  visible 
by  a  bouncing  of  the  flame  brush.  The  stabilisation  of  the  flame 
brush  off  the  symmetry  plane  of  the  burner  was  also  noticeable. 
It  was  suggested  that  an  interaction  of  the  jets  with  the  surround¬ 
ing  air  caused  both  phenomena  and  difficulties  were  reduced  with 
the  introduction  of  co-flowing  streams  as  confirmed  by  a  more 
homogeneous  turbulence  distribution  with  less  energy  at  lower 
frequencies.  Subsequent  studies  (e.g.,  [12,17])  have  shown  that 
the  flow  symmetry  can  be  much  improved  for  nozzle  separations 
H/D  ~  1,  where  H  is  the  nozzle  separation  and  D  the  nozzle 
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diameter,  with  low  frequency  instabilities  removed  and  overall  im¬ 
proved  stability  [6], 

The  comparatively  low  turbulence  levels  achieved  in  earlier 
studies  at  the  point  of  flame  extinction  have  proved  problematic 
for  computational  methods  based  on  high  Reynolds  number 
assumptions.  The  same  observation  [9]  also  applies  for  isothermal 
flows  at  Reynolds  (Re)  numbers  below  10,000  and  higher  values 
were  required  in  order  to  provide  good  agreement  with  numerical 
simulations  when  perforated  plates  were  used  as  turbulence  gen¬ 
erators.  However,  such  conditions  were  found  to  lead  to  early 
extinction  for  flames  featuring  alkane  fuels  due  to  the  presence 
of  bulk  strain  at  nozzle  separations  small  enough  (H/D  ^  1 )  to  pre¬ 
vent  axial  bulk  instabilities  in  the  flow.  Coppola  and  Gomez  [19] 
used  high  blockage  ratio  turbulence  generators  to  increase  turbu¬ 
lence  levels  substantially  and  suggested  [20]  the  use  of  the  geom¬ 
etry  as  a  benchmark  for  practical  systems.  The  turbulence 
generators  were  subsequently  found  trigger  large  scale  oscillations 
of  the  stagnation  surface  at  frequencies  <200  Hz  and  it  was  sug¬ 
gested  that  such  instabilities  need  to  be  screened,  e.g.,  using  Proper 
Orthogonal  Decomposition  (POD),  to  eliminate  artifacts  (e.g.,  in 
turbulence  levels)  from  the  measurements  [21].  However,  it  has 
subsequently  been  suggested  [22]  that  flows  with  combustion 
require  conditional  POD  (CPOD)  as  conventional  POD  does  not 
distinguish  between  fluid  structures  in  reactants  and  products. 
The  filtering  of  instabilities  using  the  latter  approach  can  hence 
arguably  only  be  applied  to  flows  without  significant  density  vari¬ 
ations.  Fractal  grids  are  here  used  to  avoid  such  difficulties  and  to 
provide  significantly  increased  turbulence  levels  without  any  is¬ 
sues  arising  with  respect  to  bulk  flow  instabilities.  The  use  of  frac¬ 
tal  grids  was  first  presented  in  the  context  of  wind  tunnels  by 
Vassilicos  and  Hunt  [23]  and  Seoud  and  Vassilicos  [24]  and  later 
by  Hurst  and  Vassilicos  [25],  The  first  application  of  fractal  grids 
in  an  opposed  jet  geometry  was  presented  by  Geipel  et  al.  [6]  in 
a  study  of  isothermal  flows.  It  was  shown  that  such  grids  can  in¬ 
crease  turbulent  Reynolds  ( Ret )  numbers  in  excess  of  a  factor  of 
two  for  the  same  bulk  velocity  as  compared  to  conventional 
perforated  plate  turbulence  generators.  It  has  also  been  shown  that 
fractal  grids  can  substantially  enhance  turbulent  diffusion  [26,27], 
The  consequence  is  a  significant  change  in  the  balance  of  turbulent 
to  bulk  strain  while  also  maintaining  acceptable  flow  symmetry 
[6[.  The  fractal  grid  approach  was  subsequently  used  by  Goh 
et  al.  [28]  to  explore  the  transition  of  premixed  JP-10  flames  from 
the  corrugated  flamelet  regime  to  a  Homogeneous  Charge  Diffu¬ 
sion  Ignition  (HCD1)  [29]  related  flameless  combustion  mode  and 
by  Goh  et  al.  [22]  to  analyse  conventional  opposed  jet  flames. 

The  current  work  extends  past  studies  of  opposed  jet  flames  by: 
(i)  The  use  of  multiscale  fractal  grid  generated  turbulence  in  a  sub¬ 
stantially  revised  configuration  aimed  at  providing  increased  tur¬ 
bulent  strain  while  providing  simplified  upstream  boundary 
conditions  in  order  to  facilitate  computational  studies.  The  turbu¬ 
lent  Reynolds  number  range  is  moved  from  50-120  to  130-318, 
as  compared  to  conventional  perforated  plate  generators,  (ii) 
Velocity  and  reaction  progress  variable  statistics,  including  condi¬ 
tional  velocities  and  scalar  fluxes,  are  reported  for  stoichiometric 
and  fuel  lean  methane,  ethylene  and  propane  flames  approaching 
extinction,  (iii)  The  turbulence  structure  obtained  using  conven¬ 
tional  and  fractal  grids  was  analysed  using  POD  and  CPOD  [22] 
techniques,  (iv)  Probability  density  functions  were  determined 
for  the  instantaneous  location  of  the  stagnation  point  and  show 
the  movement  to  be  of  the  order  of  the  integral  length  scale  and 
hence  not  influenced  by  unstable  bulk  flow  motion,  (v)  Probability 
density  functions  of  flame  curvature  were  determined  using  a 
validated  multi-step  flame  front  detection  algorithm  [22,28],  (vi) 
Estimates  of  turbulent  burning  velocities  are  provided  and  differ¬ 
ences  associated  with  alternative  determination  methods  quanti¬ 
fied.  Finally,  (vii)  the  rate  of  dissipation  was  estimated  using  the 


velocity  gradient  based  technique  of  George  and  Hussein  [30], 
Overall,  the  data  sets  present  an  opportunity  for  a  systematic  eval¬ 
uation  of  calculation  methods  for  premixed  turbulent  flames 
approaching  extinction. 

2.  Experimental  configuration,  techniques  and  uncertainties 

The  opposed  jet  geometry  used  as  a  starting  point  in  the  current 
work  consists  of  two  nozzles  in  a  vertical  arrangement  originally 
designed  by  Geyer  et  al.  [5,31]  and  is  identical  to  that  described 
by  Geipel  et  al.  [6].  Both  nozzles  were  water-cooled  to  prevent  dif¬ 
ferences  in  the  reactant  densities  due  to  preheating.  The  maximum 
positional  uncertainties  for  the  nozzles  with  the  current  configura¬ 
tion  are  0.2  mm  for  the  coaxial  alignment  and  0.5°  in  the  angular 
alignment.  The  outlet  of  each  nozzle  is  30  mm  in  diameter  and, 
for  the  base  case  configuration,  turbulence  is  generated  50  mm  up¬ 
stream  of  the  nozzle  exit  plane  using  perforated  plates  with  a  hole 
diameter  of  4  mm  and  a  blockage  of  45%  as  shown  in  Fig.  1.  The 
perforated  plates  are  similar  to  those  used  by  Mastorakos  et  al. 
[32-34],  Sardi  et  al.  [15,35,36]  and  Lindstedt  et  al.  [9,12],  Fractal 
grids  [6,23-25]  subsequently  replaced  the  perforated  plates  and 
were  used  to  increase  turbulence  intensities  at  the  nozzle  exits. 
Following  an  extensive  experimental  study  of  the  corresponding 
isothermal  flow  field  [6],  a  space-filling  fractal  cross  grid  with  a  to¬ 
tal  blockage  of  65%  and  fractal  dimensions  corresponding  to  a  max¬ 
imum  bar  width  of  2.0  mm  and  a  minimum  width  of  0.5  mm  was 
chosen  as  shown  in  Fig.  1.  The  selected  fractal  grid  increased  the 
turbulent  Reynolds  number  range  from  50-120  to  130-318  for 
the  same  bulk  velocity  range  of  4.0-8. 0  m/s,  based  on  the  integral 
length  scale  ( lt )  of  ~3.1  mm  determined  by  Geipel  et  al.  [6],  with 
turbulence  levels  measured  at  the  nozzle  exits.  Goh  [37]  deter¬ 
mined  the  corresponding  integral  length  scales  for  the  current  re¬ 
vised  geometry,  discussed  below,  and  obtained  marginally  higher 
values  of  ~3.5  mm  at  4.0  m/s  bulk  velocity  and  ~3.2  mm  at 
8.0  m/s.  The  estimation  of  length  scales  and  the  impact  on  the  fit¬ 
ting  of  longitudinal  energy  spectra  based  on  hot  wire  anemometry 
and  PIV  data  has  been  discussed  by  Geipel  et  al.  [6],  The  axial  insta¬ 
bility  of  the  flow  observed  at  higher  nozzle  separations  (e.g., 
[13,17,18])  was  effectively  removed  by  the  selected  H/D  (=1)  ratio 
and  further  reduced  by  the  presence  of  a  co-flow  [10].  The  co-flow- 
ing  stream  of  air  was  set  to  a  velocity  of  >10%  of  the  corresponding 
bulk  flow  velocity  of  the  reactant  mixture  in  order  to  reduce  the  ef¬ 
fect  of  the  shear  layer  that  forms  between  the  opposed  jet  flows 
and  the  ambient  air. 

Dry  and  filtered  air  was  supplied  to  each  nozzle  at  4  bar(g)  by  a 
compressor  using  digital  mass  flow  controllers  (Bronkhorst  UK 
LTD).  The  deviation  of  each  mass  flow  controller  was  ^0.8%  RD 
(reading)  plus  $:0.2%  FS  (full  scale)  and  the  reproducibility  better 
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Fig.  1.  Left:  Traditional  perforated  plates  at  a  blockage  ratio  of  45%.  Right: 
Dimensions  of  the  fractal  cross  grid  used  in  the  current  study.  The  dimensions 
correspond  to  Grid  I  used  in  the  study  by  Geipel  et  al.  [6]  and  provides  a  blockage 
ratio  of  65%  with  maximum  and  minimum  fractal  bar  widths  of  2.0  mm  and 
0.5  mm. 
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Fig.  2.  Opposed  jet  burner  geometry  showing  the  locations  of  the  turbulence 
generators  and  impact  plates.  The  current  geometry  is  configured  for  four  4  mm 
inlet  jets.  Alternatively,  four  8  mm  inlet  jets  were  also  used. 

than  ^0.1%.  Two  additional  mass  flow  controllers,  each  calibrated 
for  methane,  propane  and  ethylene,  were  used  to  control  the  fuel 
flows.  The  fuels  were  supplied  using  bottles  of  either  methane 
(99%  purity),  propane  (96%  purity)  and  ethylene  (99.9%  purity)  at 
3.0  bar(g).  The  mass  flow  rates  were  monitored  using  a  custom 
written  LabView  interface  connected  to  a  Readout/Control  unit 
Type  E-7000  that  was  connected  via  a  flow-bus  system  to  each  of 
the  four  mass  flow  controllers.  The  mixing  of  the  reactants  fol¬ 
lowed  the  approach  described  by  Lindstedt  et  al.  [12]  with  overall 
bulk  velocities  in  the  range  4.0-8.0  m/s. 

A  Particle  Image  Velocimetry  (PIV)  system  (LaVision  Flowmas- 
ter  3)  was  used  to  measure  the  axial  and  radial  velocity  compo¬ 
nents.  The  central  plane  perpendicular  to  the  nozzle  exits  was 
illuminated  by  two  120  mj  Solo-New  Wave  double  pulse  Nd:YAG 
lasers.  The  laser  head  was  equipped  with  LaVision  light  sheet  op¬ 
tics  that  allowed  an  adjustable  light  sheet  thickness  between 
0.5  mm  and  2.5  mm.  The  thickness  was  set  to  ~0.8  mm  for  all  mea¬ 
surements  presented.  The  flow  field  was  viewed  using  a  12  bit, 
1376  x  1040  pixel,  10  Hz  CCD  Imager  Intense  camera  equipped 
with  a  50  mm  Nikon  lens.  The  time  between  the  two  PIV-images 
was  adjusted  between  15  and  40  ps  according  to  the  mean  bulk 
velocity  in  order  to  reduce  spurious  vectors.  For  each  correlation, 
two  images  were  divided  into  smaller  interrogation  windows  and 
a  decreasing  window  size  starting  from  128  x  128,  via  64  x  64 


down  to  32  x  32  pixels  was  chosen  with  a  50%  overlap  leading  to 
a  vector  spacing  of  0.4  mm.  Host-Madsen  and  Nielsen  [38]  esti¬ 
mated  that  an  interrogation  window  size  of  0.21,  would  result  in 
approximately  a  10%  uncertainty  in  the  fluctuations.  The  current 
integral  length  scales  result  in  an  estimated  uncertainty  in  the  tur¬ 
bulent  fluctuations  of  ^15%.  All  scalar  and  velocity  data  was  ob¬ 
tained  by  averaging  at  least  1000  statistically  independent 
instantaneous  vector  fields  using  a  purpose  written  FORTRAN  pro¬ 
gram.  Further  post-processing  featured  a  window  based  algorithm 
for  the  determination  of  the  instantaneous  stagnation  point  loca¬ 
tion.  The  window  size  was  set  to  36  times  the  vector  spacing 
(0.4  mm).  The  rotation  of  the  instantaneous  stagnation  plane  was 
determined  by  applying  a  linear  fit  to  the  two  outermost  points 
of  the  stagnation  point  streamlines. 

Silicon  oil  droplets  with  an  average  diameter  of  1  pm  produced 
by  two  PALAS  Aerosol  Generators  Type  AGF  10.0  were  used  to 
determine  the  reaction  progress  variable  statistics  and  conditional 
reactant  velocities.  Velocity  data  was  also  obtained  using  alumin¬ 
ium-oxide  particles  with  an  average  diameter  of  3-5  pm  as  mea¬ 
sured  by  a  Cilas  1064  particle  size  analyser.  The  number  of 
excited  pixels  on  the  CCD-image  was  kept  between  two  and  five. 
This  is  in  keeping  with  the  optimal  particle  image  size  and  there¬ 
fore  avoided  errors  associated  with  peak-locking,  described  in 
greater  detail  by  Angele  and  Muhammad-Klingmann  [39].  Further¬ 
more,  the  RMS-uncertainty  was  kept  at  a  minimum  by  using  a 
particle  image  diameter  of  more  than  one  pixel  [40],  Other  sugges¬ 
tions  for  the  optimal  particle  image  diameter  range  from  1 .5  pixels 
[41  ],  2.0  pixels  [42]  and  up  to  4.0  pixels  [43 ].  The  latter  value  was 
obtained  by  a  Monte-Carlo  calculation  including  velocity  gradients 
due  to  shear.  The  authors  suggested  that  the  presence  of  shear  was 
the  reason  for  the  increase  of  the  optimal  particle  size  on  the  PIV- 
images. 

Stagnating  flows  have  a  high  dynamic  range  with  the  minimum 
velocity  approaching  zero  at  the  stagnation  plane  where  the  largest 
errors  can  be  expected  [12],  The  relaxation  time  of  a  seeding  par¬ 
ticle  to  a  step  change  in  the  fluid  velocity  was  estimated  using 
the  approach  of  Han  and  Mungal  [44], 


(1) 


The  relaxation  time  for  the  dp  =  1  pm  silicon  oil  droplets 
(pp  =  980  kg/m3)  was  estimated  at  ~3  ps  and  for  the  dp  =  3  pm 
A1203  particles  ( pp  =  3900  kg/m3)  the  relaxation  time  was  ~29  ps 
in  the  products  (p  =  6.75941CT5  kg/m/s)  and  ~109  ps  in  the  reac¬ 
tants  (p  =  1.78351CT5  kg/m/s)  [45],  The  smallest  PIV  scale  based 
on  the  largest  bulk  velocity  of  8.0  m/s  and  the  smallest  interroga¬ 
tion  window  size  (32  x  32  pixels)  is  around  50  ps  and  sufficiently 
large  for  particles  to  follow  the  flow  in  the  products.  Particle  lag  in 
the  reactants  could  impact  the  A1203  seeding  though  previous 
work  [12]  has  shown  that  the  reactant  velocities  obtained  for  both 
types  of  seeding  were  within  15%. 

Optical  measurement  techniques  that  rely  on  seeding  particles 
(e.g„  LDV  and  PIV)  can  suffer  from  effects  of  thermophoresis  lead¬ 
ing  to  a  particle  velocity  lag  as  shown  by  Sung  et  al.  [46,47].  Gomez 
and  Rosner  [48]  suggest  a  local  particle  drift  velocity  (vT)  as  shown 
Eq.  (2),  which  was  approximated  by  Amantini  et  al.  [49]  using  Eq. 
(3),  where  aT  is  the  dimensionless  thermophoretic  diffusion  factor, 
D  the  particle  Brownian  diffusivity,  T  the  temperature,  am  the  tan¬ 
gential  momentum  accommodation  coefficient  between  the  fluid 
and  the  solid  particle  (often  assumed  to  be  unity)  and  v  the  kine¬ 
matic  viscosity. 


-VT 

vT  =  (<xT  D)partjcle  -jr- 
j/T~0.75[l+|am]v 


-VT 


(2) 

(3) 
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In  the  vicinity  of  the  flame  front  the  effect  of  thermophoresis 
might  be  greater  than  the  uncertainty  in  the  PIV  data  in  which  case 
the  effect  needs  to  be  considered  as  shown  by  Amantini  et  al.  [49], 
Chelliah  et  al.  [50]  and  Sung  et  al.  [46],  For  an  ethylene  flame 
(</>  =  0.9)  it  was  estimated  that  vT  ~  0.49  m/s  based  on  premixed 
laminar  flame  data.  The  resulting  uncertainty  in  velocity  statistics 
across  the  flame  front  is  ^20%  when  evaluated  relative  to  the  corre¬ 
sponding  bulk  velocity.  For  other  flames  the  uncertainty  was  ^10%. 

Beam  steering  effects  caused  by  a  change  in  the  refractive  index 
of  the  fluid  and  the  impact  of  the  density  gradient  across  a  flame 
front  were  quantified  by  Kaiser  et  al.  [51  j  using  a  ray-tracing  tech¬ 
nique  based  on  Rayleigh  images  as  an  approximate  measure  of  the 
index-of-refraction  field.  The  study  indicated  that  the  spread  of  the 
rays  in  the  image  plane  was  on  the  sub-pixel  scale  and  therefore 
negligible  as  also  discussed  by  Frank  et  al.  [52],  Bohm  [53]  found, 
by  comparing  time-correlated  OH-LIF  images  obtained  in  an  op¬ 
posed  jet  geometry,  that  the  flame  movement  between  two  images 
50  ps  apart  was  negligible  and  hence  the  uncertainty  of  the  veloc¬ 
ity  measurements  caused  by  beam  steering  in  the  measured  time 
interval  could  be  ignored.  The  time  interval  between  the  PIV 
images  in  the  current  work  was  between  15  and  40  ps  and  uncer¬ 
tainties  caused  by  beam  steering  were  neglected. 

2.1.  Flame  front  detection 

The  reaction  progress  variable  statistics  were  produced  from 
PIV  images  using  a  custom  written  density  segregation  algorithm, 
described  by  Goh  et  al.  [22,28],  to  identify  gradients  created  by 
changes  in  seeding  density  due  to  heat  release.  Similar  algorithms 
have  been  developed  and  shown  capable  of  producing  excellent 
agreement  with  flame  fronts  obtained  using  simultaneous  CH- 
LIF/PIV  by  Steinberg  et  al.  [54]  or  0H/CH20-L1F/PIV  by  Pfadler 
et  al.  [55].  The  current  algorithm  was  initially  validated  [22]  using 
images  of  a  silicon  oil  seeding  that  evaporates  as  it  passes  through 
the  flame  front  as  previously  discussed  for  stagnating  flows  by 
Shepherd  et  al.  [57],  Kostiuk  et  al.  [13]  and  Lindstedt  et  al.  [12], 
The  technique  was  further  validated  by  Goh  et  al.  [28]  in  the  pres¬ 
ence  of  reaction  zone  broadening  by  comparing  with  OH-PLIF 
images  obtained  using  the  Ba(N03)2  crystal  based  technique  of  Kerl 
et  al.  [56],  The  first  step  in  the  algorithm  is  to  determine  the  mean 
intensity  of  the  Mie-scattering  in  a  defined  reference  window.  Goh 
et  al.  [22]  have  discussed  the  thresholding  process  in  detail  and 
inspection  of  particle  images  suggest  a  Gaussian  nature.  The  step 
naturally  filters  out  background  noise  should  images  contain  sig¬ 
nificant  noise  from  reflections,  non-uniform  laser  intensity  or  laser 
sheet  thickness.  The  filtering  produces  binary  images  with  an  equal 
weight  given  to  all  accepted  peaks.  Subsequently,  using  the  density 
of  peaks  in  the  reference  region,  gaps  between  peaks  were  filled 
using  a  moving  average  filter  to  avoid  leaving  gaps  in  the  reactant 
regions  due  to  minor  fluctuations  in  seeding  densities.  Subse¬ 
quently,  a  bilateral  filter  [58]  was  used  to  reduce  the  intensity  fluc¬ 
tuations  in  the  reactant  and  product  regions  while  preserving  the 
steep  gradients  at  the  flame  front.  The  latter  was  achieved  by  using 
convolution  of  range  and  distance  filters  with  the  image  produced 
by  the  moving  average  filter.  The  resulting  image  had  a  relatively 
smooth  histogram  with  multiple  peaks  representing  segmented  re¬ 
gions  of  similar  intensities.  By  applying  a  cubic  spline  fit  to  the  his¬ 
togram  of  the  resulting  image,  as  applied  by  Pfadler  et  al.  [59]  to 
filtered  images,  the  nearest  minima  of  the  intensity  lower  than 
the  minimum  intensity  in  the  reference  window  was  used  as  a 
threshold  for  creating  a  binary  image.  Canny  edge  detection  [60] 
with  a  Sobel  operator  was  then  applied  on  the  resulting  binary  im¬ 
age  to  obtain  steep  gradients.  Subsequently,  non-maxima  suppres¬ 
sion,  followed  by  a  skeleton  generating  algorithm  [61],  were  used 
to  thin  the  resulting  edges  to  a  thickness  of  one  pixel  with  the 
flame  front  identified  as  the  longest  line  closest  to  the  reference 


window.  The  process  avoided  the  identification  of  the  flame  front 
from  the  opposing  reactant  stream  should  seeding  densities  in 
both  reference  windows  be  the  same. 

Two  data  sets,  obtained  with  A1203  or  silicon  oil  seeding  with 
each  containing  1000  PIV  images,  were  obtained  for  identical  flow 
conditions.  Examples  of  images  are  shown  in  Fig.  3  where  each  pix¬ 
el  is  around  0.025  mm  in  size.  Silicon  oil  evaporates  across  the 
instantaneous  flame  front  resulting  in  a  binary  PIV  image  without 
further  treatment.  The  developed  algorithm  was  applied  to  both 
sets  of  data  and  the  first  two  moments  of  the  reaction  progress  var¬ 
iable  statistics  showed  only  subtle  differences  as  exemplified  in 
Fig.  4  for  methane/air  flames  at  </>  =  0.90.  The  slight  asymmetry  ob¬ 
served  is  of  the  order  0.4  mm  and  close  to  the  estimated  Kolmogo¬ 
rov  scale  in  the  combustion  products.  The  mean  progress  variable 
was  marginally  lower  for  aluminium  oxide  in  the  vicinity  of  the 
stagnation  point.  However,  the  effect  is  arguably  due  to  the  greater 
heat  losses  caused  by  the  aluminium  oxide  seeding  which  results 
in  lower  flame  temperatures  and  a  modest  decrease  in  the  distance 
between  the  twin  flames.  The  current  silicon  oil  evaporates  at 
550  K  and  the  detected  flame  fronts  will  be  biased  towards  the 
reactant  streams.  The  discrepancies  are  of  the  order  of  0.5  mm, 
which  is  comparable  to  the  laminar  flame  thickness  as  determined 
by  simulations  featuring  detailed  chemistry  [45].  Flame  curvature 
statistics  were  also  calculated  to  provide  a  further  example.  The 
x(s)  and  y(s)  coordinates,  where  s  denotes  the  distance  along  the 
curve,  were  fitted  as  close  to  the  normalised  curve  as  possible 
using  cubic  splines.  Subsequently,  the  curvature  h(s )  was  calcu¬ 
lated  using  Eq.  (4)  [64], 


h(s) 


dx  d2y d2x  dy 

ds  ds7  ds7  ds 


(4) 


Based  on  the  coordinate  system  the  curvature  is  positive  for 
upper  flame  fronts  convex  to  products  and  for  lower  flame  fronts 
concave  to  products.  The  use  of  cubic  splines  was  found  necessary 
as  traditional  (e.g.,  finite  difference)  methods  created  noise  in  the 
curvature  calculations.  The  curvature  statistics  presented  here 
were  solely  used  for  verification  of  the  algorithm  as  precise  curva¬ 
ture  calculations  across  a  large  range  of  lengthscales  remains 
intractable.  The  calculated  statistics  from  500  images  in  each  set 
are  presented  in  Fig.  5.  There  is  agreement  for  the  curvature  pdf 
for  the  upper  flame,  with  some  discrepancies  for  the  lower  flame 
close  to  zero  curvature.  The  data  presented  were  obtained  with  a 
perforated  plate  without  a  fitted  impact  plate  and  individual 
frames  were  investigated  to  ensure  that  the  effect  was  not  caused 
by  the  flame  front  detection  algorithm. 


2.2.  Conditional  proper  orthogonal  decomposition 


The  instantaneous  progress  variable  (c)  was  also  used  to  segre¬ 
gate  the  reactants  from  the  products  to  obtain  unconditional  ( u ) 
Proper  Orthogonal  Decomposition  (POD)  and  conditional  (CPOD) 
reactant  (r)  and  product  (p)  data  using  the  method  of  Goh  et  al. 
[22]  with  the  spatial  mean  of  the  turbulence  kinetic  energy  (£)  ob¬ 
tained  using  Eqs.  (5). 
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Fig.  3.  Verification  of  the  flame  front  detection  algorithm  using  aluminium  oxide  seeding  (left)  and  silicon  oil  droplets  (right).  Lines  indicate  detected  flame  fronts. 
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Fig.  4.  A  comparison  of  centreline  results  for  the  first  and  second  moments  of  the 
reaction  progress  variable  (c)  obtained  using  the  flame  front  detection  algorithm  in 
combination  with  aluminium  oxide  and  silicon  oil  seeding.  Solid  line:  Aluminium 
oxide  seeding.  Dashed  line:  Silicon  oil  seeding. 
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Fig.  5.  The  pdf  of  curvature.  Bars:  Aluminium  oxide  seeding.  Lines:  Silicon  oil 
seeding.  Left:  Upper  flame  front.  Right:  Lower  flame  front. 


The  spatial  mean  of  the  absolute  kinetic  energy  represented  by 
each  mode  was  determined  by  simply  calculating  £,  the  normalised 
values  £m  n0rm  and  the  cumulative  energy  represented  by  M  modes. 

Goh  et  al.  [22]  also  applied  the  method  of  George  and  Hussein 
[30]  to  estimate  the  rate  of  turbulence  dissipation  (e)  from  PIV  data 
using  a  linear  combination  of  mean  square  values  of  local  velocity 
gradients  and  the  kinematic  velocity  (v).  The  Kolmogorov  length 
scale  was  estimated  to  be  ~0.1  mm  in  the  reactants  and 
~0.5  mm  in  the  products  under  similar  conditions.  The  corre¬ 
sponding  Kolmogorov  velocities  were  estimated  to  be 


~0.16  ms  1  in  the  reactants  and  ~0.78  ms  1  in  the  products.  The 
same  method  was  applied  in  the  current  work  as  discussed  below. 

3.  Results 

The  current  study  presents  data  obtained  for  turbulent  pre¬ 
mixed  opposed  jet  flames  featuring  methane,  propane  and  ethyl¬ 
ene  fuels  as  well  as  isothermal  flow  conditions.  The  latter  were 
included  to  provide  additional  information  on  the  flow  at  the  noz¬ 
zle  exits  as  well  as  further  information  on  the  influence  of  the  flow 
path  upstream  of  the  turbulence  generators  on  velocity  statistics 
and  eigenmodes  as  determined  using  POD  and  CPOD. 

The  section  upstream  of  the  turbulence  generators  was  modi¬ 
fied  by  fitting  four  4  mm  or  8  mm  equally  spaced  inlet  jets  to  ex¬ 
plore  the  impact  of  flow  conditions  leading  up  to  the  turbulence 
generators.  For  the  same  purpose,  additional  impact  plates  were 
also  placed  upstream  of  the  turbulence  generating  plates  in  order 
to  break  up  the  four  inlet  jets  and  to  explore  the  impact  on  the 
resulting  turbulence  intensities  at  the  nozzle  exits.  The  burner 
geometry  is  shown  in  Fig.  2  and  the  six  investigated  turbulence 
generator  configurations  are  summarised  below. 

•  Case  1:  Conventional  turbulence  generator  without  fitted 
impact  plate  and  (a)  4  mm  or  (b)  8  mm  equally  spaced  inlet  jets. 

•  Case  2:  Conventional  turbulence  generator  with  fitted  impact 
plate  and  (a)  4  mm  or  (b)  8  mm  equally  spaced  inlet  jets. 

•  Case  3:  Fractal  grid  without  fitted  impact  plate  and  4  mm 
equally  spaced  inlet  jets. 

•  Case  4:  Fractal  grid  with  fitted  impact  plate  and  4  mm  equally 
spaced  inlet  jets. 

The  above  configurations  were  initially  investigated  for  isother¬ 
mal  flows  in  preparation  for  the  subsequent  studies  featuring 
combustion. 

3.1.  Isothermal  flows 

The  mean  axial  and  radial  velocity  profiles  obtained  using  tradi¬ 
tional  perforated  plate  turbulence  generators  (Cases  1  and  2)  under 
isothermal  flow  conditions  are  presented  in  Fig.  6  for  a  location 
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1.5  mm  downstream  of  the  nozzle  exit  plane.  The  current  mean 
velocity  profiles  show  the  same  trends  as  presented  by  Lindstedt 
et  al.  [9 1  and  Geyer  et  al.  [31  ].  A  reduction  is  apparent  for  all  mean 
axial  velocity  profiles  towards  the  burner  centreline  due  to  the 
pressure  increase  along  the  stagnation  point  streamline,  as  also 
discussed  by  Korusoy  and  Whitelaw  [8],  for  all  four  upstream  flow 
conditions.  The  use  of  the  additional  impact  plates  (Cases  2a  and 
2b)  resulted  in  a  stronger  velocity  gradient  towards  the  nozzle  cen¬ 
tre  for  the  4  mm,  as  compared  to  the  8  mm,  jets.  Further  compar¬ 
isons  of  the  nozzle  exit  profiles  show  that  the  impact  plates  also 
suppressed  the  influence  of  the  nozzle  inlet  velocity  and  ensured 
independence  of  the  mean  velocity  exit  profiles  of  conditions  fur¬ 
ther  upstream.  However,  axial  and  radial  velocity  fluctuations 
across  the  nozzle  exits  exhibit  a  much  greater  dependency  on  the 
flow  path  upstream  of  the  turbulence  generators  as  shown  in 
Fig.  7.  The  use  of  4  mm  inlet  jets  without  impact  plates  provided 
turbulence  fluctuations  up  to  20%  of  the  corresponding  bulk  veloc¬ 
ity  (U b)  as  compared  to  values  around  9%  with  impact  plates  fitted. 
The  latter  value  is  in  agreement  with  data  presented  in  earlier 
studies  [6,9  .  An  integral  length  scale  of  2.6  mm  for  the  perforated 
plates  [12]  results  in  a  turbulent  Reynolds  number  of  approxi¬ 
mately  50  for  the  flow  using  impact  plates.  The  increased  velocity 
fluctuations  caused  by  the  removal  of  the  impact  plates  is  solely 
the  result  of  direct  jet  impingement  on  the  perforated  plates.  A 
reduction  of  the  inlet  jet  velocity  by  a  factor  of  four,  through  the 
use  of  8  mm  jets,  resulted  in  a  100%  decrease  in  the  velocity  fluctu¬ 
ations  across  the  nozzle  diameter.  Accordingly,  it  can  be  concluded 
that  the  conventional  perforated  plate  turbulence  generators  are 
strongly  dependent  upon  the  upstream  flow  conditions. 

Exchanging  the  perforated  plates  with  fractal  grids  provided  the 
expected  ~100%  increase  in  fluctuation  levels  [6].  The  velocity  sta¬ 
tistics  1 .5  mm  downstream  of  the  nozzle  exits  are  presented  in 
Fig.  8  for  Cases  3  and  4.  Due  to  (small)  uncertainties  in  the  manu¬ 
facturing  process  of  the  fractal  grid  structure,  it  cannot  be  assumed 
a  priori  that  the  nozzle  exit  profiles  are  completely  symmetric  in 
the  circumferential  direction.  Therefore,  velocity  statistics  for  each 
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Fig.  6.  Axial  and  radial  mean  velocity  profiles  plotted  across  the  nozzle  radius 
1.5  mm  downstream  of  the  nozzle  exit.  Upper  row:  4  mm  jets  and  a  bulk  velocity  of 
4.0  m/s.  Lower  row:  8  mm  jets  and  a  average  bulk  velocity  of  4.0  m/s.  Symbols:  o 
indicate  the  use  of  impact  plates  and  a  no  impact  plates,  used  in  combination  with 
traditional  perforated  plates. 
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Fig.  7.  Influence  of  the  upstream  boundary  conditions  on  the  velocity  fluctuations 
at  a  mean  bulk  velocity  of  4.0  m/s.  Symbols:  a  indicates  4  mm  jets  and  □  8  mm  jets 
measured  without  impact  plates;  o  indicates  4  mm  jets  and  0  8  mm  jets  with  the 
use  of  additional  impact  plates,  used  in  combination  with  traditional  perforated 
plates. 


Fig.  8.  Top  row:  axial  (left)  and  radial  (right)  mean  velocity  profiles  measured 
across  the  nozzle  radius  at  a  location  1.5  mm  downstream  of  the  nozzle  exit  for  a 
bulk  velocity  of  4.0  m/s.  Lower  row:  normalised  axial  (left)  and  radial  (right) 
fluctuations;  o  indicates  the  use  of  impact  plates  and  a  without  impact  plates.  In 
both  cases  four  4  mm  inlet  jets  were  used  in  combination  with  fractal  grids. 


flow  condition  were  obtained  using  four  separate  measurements 
each  with  a  45°  rotation  in  the  circumferential  direction  and  the 
error  bars  shown  in  each  graph  indicate  the  maximum  deviations 
obtained  for  the  four  data  sets.  Differently  from  the  conventional 
grids,  it  is  apparent  that  the  introduction  of  impact  plates  does 
not  cause  a  significant  reduction  in  turbulence  intensities.  Compar¬ 
isons  of  the  velocity  statistics  obtained  using  4  mm  inlet  jets  with 
and  without  the  additional  impact  plates  combined  with  the  con¬ 
ventional  turbulence  generator  (Cases  la  and  2a)  and  the  fractal 
grid  (Cases  3  and  4)  are  shown  in  Fig.  9  for  the  mean  axial  velocities 
and  in  Fig.  10  for  the  axial  and  radial  fluctuations.  A  reduced  mean 
axial  velocity  was  observed  for  both  turbulence  generators  when 
additional  impact  plates  were  used.  However,  no  significant  de¬ 
crease  in  turbulence  fluctuations  were  found  when  combining 
fractal  based  turbulence  generation  with  impact  plates.  Accord¬ 
ingly,  the  fractal  grids  have  the  notable  advantage  of  removing 
the  sensitivity  to  upstream  conditions  while  maintaining  signifi¬ 
cantly  increased  levels  of  turbulence. 

The  flow  field  differences  outlined  above  were  further  analysed 
by  applying  a  POD  technique  within  the  region  r  =  ±D/2  to 
decompose  the  instantaneous  vector  fields  into  their  empirical 
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Furthermore,  a  mode  switch,  as  compared  to  the  conventional  tur¬ 
bulence  generators,  was  observed.  The  lower  proportion  of  the 
kinetic  energy  contained  in  the  first  and  second  modes  for  the  frac¬ 
tal  grid  turbulence  generator  suggests  that  the  turbulence  intensi¬ 
ties  of  «20%  of  the  bulk  velocity  at  the  nozzle  exit  plane  are  not 
dominated  by  the  first  two  modes.  The  results  for  fractal  grids 
are  hence  unlike  traditional  perforated  plates  at  comparable  levels 
of  fluctuations  and  generate  turbulence  with  a  different  structure 
by  exciting  a  larger  number  of  modes.  The  impact  upon  premixed 
turbulent  flames  is  explored  below. 


Fig.  9.  Mean  normalised  axial  velocity  profiles  along  the  burner  centreline  for  a 
bulk  velocity  of  4.0  m/s  and  using  4  mm  inlet  jets  with  zero  indicating  the 
stagnation  point.  Symbols:  a  indicates  results  obtained  without  impact  plates  and  o 
with  impact  plates.  Left:  Traditional  perforated  plates;  Right:  Fractal  grids. 


Fig.  10.  Normalised  axial  and  radial  velocity  fluctuations  along  the  burner 
centreline  for  a  bulk  velocity  of  4.0  m/s  and  with  4  mm  inlet  jets  with  zero 
indicating  the  stagnation  point.  Symbols:  a  indicates  results  obtained  without 
impact  plates  and  o  with  impact  plates.  Left:  Traditional  perforated  plates;  Right: 
Fractal  grids. 


eigenfields  in  order  to  explore  the  flow  dynamics.  The  method  is 
explained  in  detail  by  Holmes  et  al.  [62]  and  was  successfully  ap¬ 
plied  amongst  others  by  Adrian  et  al.  [63],  who  showed  the  effec¬ 
tiveness  of  POD  as  a  filtering  tool  using  P1V  measurements.  The 
first  eigenmode  obtained  with  the  traditional  perforated  plates 
without  (Case  la)  and  with  (Case  2a)  impact  plates  shows  a  rota¬ 
tional  motion  close  to  the  stagnation  point  that  causes  the  stagna¬ 
tion  plane  to  tilt.  The  second  mode  is  an  axial  movement  causing  a 
stagnation  plane  shift,  as  previously  analysed  by  Geipel  et  al.  [6], 
The  first  five  eigenmodes  contain  60%  of  the  total  kinetic  energy 
without  impact  plates  and  39%  with  impact  plates.  The  kinetic  en¬ 
ergy  in  the  first  mode  is  reduced  from  31.2%  to  17.6%  and  second 
mode  from  15.5%  to  11.1%  with  fitted  impact  plates.  Hence,  the  in¬ 
crease  in  the  fluctuations  using  traditional  perforated  plates  is 
more  likely  to  be  caused  by  coherent  flow  structures  generated 
by  the  four  high  velocity  inlet  jets  impinging  on  the  turbulence 
generator. 

The  fractal  generated  turbulence  without  impact  plates  (Case  3) 
contained  20.2%  of  the  kinetic  energy  in  the  first  mode  and  16.8%  in 
the  second  mode.  The  additional  impact  plates  (Case  4)  reduced 
the  first  mode  to  about  17.5%  and  the  second  to  about  12.8%. 


3.2.  Combusting  flows 

Velocity  and  scalar  statistics  are  presented  with  a  focus  on  the 
burner  axis  and  the  symmetry  plane  between  the  nozzles.  The 
eight  lean  methane,  propane  and  ethylene  flames  cover  a  range 
of  stoichiometries  (0.70  ^  <f>  ^  1 .0)  and  extend  previous  work  of 
Lindstedt  et  al.  [12]  to  cover  increased  turbulent  intensities  while 
maintaining  the  bulk  strain  rates.  Results  from  additional  vector 
field  analyses  and  the  flame  front  detection  algorithm  are  also  pre¬ 
sented  alongside  conditional  and  unconditional  velocity  statistics, 
with  Q  indicating  a  conventional  and  Q  a  Favre  average  for  the 
property  Q,  and  an  evaluation  of  alternative  methods  for  the  deter¬ 
mination  of  the  turbulent  burning  velocity.  The  structure  of  the 
reacting  turbulent  flow  field  is  further  analysed  using  the  POD 
and  CPOD  techniques  outlined  above. 

3.2.1.  Analysis  of  the  flow  field 

The  2D  detection  algorithm  developed  by  Geipel  et  al.  [6]  was 
used  to  determine  the  instantaneous  stagnation  point  location 
and  to  quantify  the  flame  movement  caused  by  turbulence  and 
any  low  frequency  bulk  flow  instabilities.  Three  examples  of  stag¬ 
nation  point  location  statistics  and  the  pdf  for  the  overall  tilt  of  the 
stagnation  plane  are  given  in  Fig.  11  for  premixed  methane/air 
flames  with  </>  =  0.9.  The  stagnation  point  essentially  moves  in  an 
envelope  of  4  x  5  mm  in  the  axial  and  radial  directions.  The  results 
suggest  that  the  turbulence  intensity  (cf.  Fig.  10)  impacts  the  width 
of  the  pdf  with  a  substantial  reduction  observed  for  the  perforated 
plates  with  a  fitted  impact  plate.  The  fractal  grid  provides  a  nar¬ 
rower  pdf  than  the  perforated  plate  for  similar  turbulence  levels. 
The  axial  movement  is  in  good  agreement  with  findings  for  the  iso¬ 
thermal  flow  field  [6],  while  the  movement  in  the  radial  direction 
increased  by  up  to  2  mm  (>50%)  due  to  combustion.  A  slight  bias 
towards  a  positive  stagnation  point  tilt  was  also  observed  for  all 
three  flow  configurations.  The  effect  was  investigated  by  compar¬ 
ing  four  different  measurement  planes  resulting  from  successive 
45°  rotations  for  the  configuration  using  fractal  grids  with  addi¬ 
tional  impact  plates  shown  in  Fig.  12.  The  stagnation  point  move¬ 
ment  for  all  four  positions  was  within  a  3  x  4  mm  envelope, 
compared  to  the  estimated  integral  length  scale  of  3.5  mm  in  the 
reactants,  suggesting  a  modest  increase  in  bias  for  the  first  mea¬ 
surement  plane  featuring  a  45°  shift  with  respect  to  the  largest 
fractal  cross. 

The  increased  radial  movement  of  the  instantaneous  stagnation 
point  was  investigated  further  by  an  analysis  of  the  first  three 
eigenmodes  for  premixed  methane  flames  with  </>  =  0.9  as  shown 
in  Fig.  13.  It  is  apparent  that  combustion  typically  modifies  or 
changes  the  dominant  eigenmodes  as  compared  to  the  isothermal 
flow.  The  first  eigenmode  for  the  case  with  conventional  perforated 
discs  without  impact  plates  contains  20.6%  of  the  total  kinetic  en¬ 
ergy  and  remains  a  rotational  motion  close  to  the  nominal  stagna¬ 
tion  point.  The  latter  is  axially  skewed  resulting  in  the  broad  pdf  of 
the  stagnation  plane  angle  shown  in  Fig.  11.  The  same  first  mode 
was  also  found  for  the  fractal  grids  with  impact  plates.  However, 
the  kinetic  energy  was  lower  at  approximately  10%  of  the  total.  A 
direct  consequence  is  a  much  narrower  pdf  of  the  rotation  angle 
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Fig.  11.  Left:  Rotation  of  the  instantaneous  stagnation  plane.  Right:  Probability  density  functions  (pdfs  in  mm  2)  of  the  instantaneous  stagnation  point  location.  Methane 
flames  with  </>  =  0.9  and  a  bulk  velocity  of  4.0  m/s.  Top:  Perforated  plates  without  impact  plates.  Middle:  perforated  plates  with  impact  plates.  Bottom:  Fractal  grids  with 
impact  plates. 


at  comparable  levels  of  fluctuations.  In  the  presence  of  combustion, 
a  new  radial  eigenmode  was  also  observed  as  shown  in  Fig.  13.  The 
associated  kinetic  energy  was  found  to  have  a  direct  influence  on 
the  fluctuation  levels  in  the  vicinity  of  the  stagnation  plane  and 
to  lead  to  an  increase  of  up  to  100%  as  compared  the  corresponding 
isothermal  flow  as  shown  in  Fig.  14.  It  has  been  shown  by  Goh  et  al. 
[22]  that  conditional  POD  (CPOD)  is  required  to  determine  the  en¬ 
ergy  distribution  and  the  hybrid  structures  obtained  using  conven¬ 
tional  POD  are  only  presented  to  show  the  relative  effects  of 
different  flow  configurations. 

Kostiuk  et  al.  [11  ]  applied  different  limiting  values  on  the  max¬ 
imum  axial  movement  of  the  stagnation  point  away  from  the  nom¬ 
inal  location  in  an  effort  to  assess  the  potential  impact  of  low 
frequency  instabilities  on  velocity  statistics  and  flow  asymmetries. 
In  the  current  work,  the  filtering  algorithm  of  Geipel  et  al.  [6]  was 
used  to  investigate  the  impact  of  stagnation  point  movement  on 
the  velocity  statistics.  The  current  filter  is  based  on  the  2D  detec¬ 
tion  algorithm  discussed  above  and  rejects  all  PIV  derived  vector 
fields  where  the  stagnation  point  location  falls  outside  a  chosen 
circular  domain  of  fixed  radius.  Filtered  velocity  statistics  for  meth¬ 
ane/air  flames  using  a  threshold  of  3.5  mm,  close  to  the  integral 
length  scale,  show  minimal  impact  on  the  axial  mean  and  fluctuat¬ 
ing  velocity  statistics  and  only  a  slight  improvement  in  symmetry. 
The  radial  fluctuations  were  reduced  by  ~7.4%  for  </>  =  0.9  and 
~4.8%  for  </>  =  0.8  at  the  nominal  stagnation  point.  A  further  reduc¬ 


tion  in  the  radial  fluctuations  was  found  when  the  threshold  was 
decreased  to  2.5  mm  though,  as  pointed  out  by  Geipel  et  al.  [6], 
it  is  probably  not  reasonable  to  apply  a  filter  less  than  the  integral 
length  scale.  Filtering  at  3.5  mm  leads  to  a  rejection  between  9.0% 
and  10.9%  of  the  PIV  images,  whereas  a  threshold  of  2.5  mm  rejects 
up  to  33.1%.  The  resulting  profiles  are  shown  in  Fig.  15  and  the 
same  trends  were  found  for  all  fuels  and  stoichiometries.  Given 
the  modest  impact  of  filtering  on  velocity  statistics,  subsequent 
graphs  present  unfiltered  data  and  the  error  bars  include  both 
measurement  uncertainties  and  flow  asymmetries. 

3.2.2.  Velocity  statistics  in  fractal  grid  generated  turbulence 

The  discussion  of  velocity  statistics  for  combusting  flows  is 
focussed  on  fractal  grid  generated  turbulence  with  fitted  impact 
plates  (Case  4)  with  selected  further  comparisons  to  highlight  dif¬ 
ferences  compared  to  perforated  plates.  Isothermal  results  show 
that  this  configuration  effectively  removes  low  frequency  flame 
movement  and  suppresses  the  impact  of  upstream  conditions 
while  producing  velocity  fluctuations  ~20%  (130  ^  Ret  ^  318)  of 
the  corresponding  bulk  velocity. 

Comparisons  of  normalised  mean  axial  velocities  along  the  bur¬ 
ner  centreline  and  corresponding  radial  velocities  along  the  stag¬ 
nation  plane  normal  to  the  nozzles  are  shown  in  Fig.  16  for 
methane  and  propane  flames  at  Ub  =  4.0  m/s  and  for  ethylene 
flames  at  Ub  =  8.0  m/s.  The  flames  cover  a  range  of  stoichiometries. 
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Fig.  12.  Left:  Rotation  of  the  instantaneous  stagnation  plane.  Right:  Probability  density  functions  (pdfs  in  mirr2)  of  the  instantaneous  stagnation  point  location.  Methane 
flames  with  o  =  0.9  and  a  bulk  velocity  of  4.0  m/s  measured  at  four  circumferential  positions.  First  and  third  row  are  measured  at  45°  and  135°  rotation  with  respect  to  the 
largest  fractal  cross.  The  second  and  last  rows  are  measured  along  the  two  perpendicular  parts  of  the  largest  fractal  cross  shown  in  Fig.  1. 


An  acceleration  in  the  axial  direction,  due  to  the  reduction  in  den¬ 
sity,  can  be  seen  for  all  fuels.  The  velocity  increase  is  less  noticeable 
for  leaner  flames  due  to  the  lower  heat  release.  An  acceleration  in 
the  radial  direction  from  zero  at  the  stagnation  point  up  to  1.25- 
1.75 Ub  at  approximately  half  a  nozzle  diameter  was  found  for  all 
fuels.  A  decrease  in  the  acceleration  of  the  radial  velocities  was  also 
evident  at  lower  stoichiometries.  A  comparison  of  radial  velocities 
with  the  corresponding  isothermal  flow  [6]  shows  that  the  mean 
radial  velocity  is  increased  by  a  factor  of  two  due  to  volume  expan¬ 
sion  caused  by  combustion. 

Axial  and  radial  fluctuations  along  the  burner  centreline  are 
shown  in  Fig.  17  for  methane  and  propane  flames  with 


Ub  =  4.0  m/s  and  for  ethylene  flames  with  Ub  =  8.0  m/s.  A  compari¬ 
son  of  the  three  fuels  at  a  stoichiometry  of  0.80  shows  that  the  nor¬ 
malised  axial  fluctuations  increase  to  a  maximum  of  0.25l/t  at  the 
stagnation  point  for  ethylene  flames,  whereas  methane  and  pro¬ 
pane  flames  showed  slightly  lower  values  of  approximately 
0.22 Ub.  All  axial  fluctuations  increase  ahead  of  the  flame  front 
and,  as  the  flow  accelerates  in  the  axial  direction,  the  turbulent 
intensities  reduce  before  increasing  again  to  reach  their  maximum 
at  the  stagnation  point.  The  complex  pattern  of  flame  and  flow  in¬ 
duced  turbulence  generation/destruction  is  characteristic  of  op¬ 
posed  jet  flames.  Radial  fluctuations  stay  constant  up  to  a 
distance  of  0.1 8D  and  increase  to  a  maximum  value  of  0.30Ub  for 
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Fig.  13.  Comparison  of  the  first  three  eigenmodes  for  premixed  methane  flames  with  </>  =  0.9  and  a  bulk  velocity  of  4.0  m/s.  Left  column:  Perforated  plates  with  no  impact 
plates.  First  eigenmode  with  20.6%,  second  eigenmode  with  14.2%  and  third  eigenmode  with  9.0%  of  the  turbulent  kinetic  energy.  Middle  column:  Perforated  plates  with 
additional  impact  plates.  First  eigenmode  with  15.7%,  second  eigenmode  with  9.3%  and  third  eigenmode  with  8.0%  of  the  kinetic  energy.  Right  column:  Fractal  grids  with 
additional  impact  plates.  First  eigenmode  with  10.0%,  second  eigenmode  with  8.9%  and  third  eigenmode  with  7.6%  of  the  kinetic  energy. 


methane  and  propane  flames  at  <j>  =  0.90.  Leaner  flames  at  a  stoichi¬ 
ometry  of  tj>  =  0.80  show  a  reduction  of  the  radial  fluctuations  by 
2.5%.  Values  between  0.26  and  0.27514  were  obtained  for  all  three 
fuels  at  </>  =  0.80,  which  suggests  that  the  radial  fluctuations  are  not 
dependent  on  the  fuel  type.  The  normalised  fluctuations  obtained 
for  the  isothermal  flows  at  the  nominal  stagnation  point  were  of 
the  order  of  0.2814  in  the  axial  and  0.1 6L4  in  the  radial  directions 
(see  Fig.  14),  respectively.  The  corresponding  values  for  the  com¬ 
busting  flow  are  0.185-0.2714  in  the  axial  and  0.24-0.30L4  radial 
directions  for  equivalence  ratios  from  0.70  to  0.90.  The  results  sug¬ 
gest  that  for  the  current  conditions,  the  axial  stresses  are  compara¬ 
ble  to  the  isothermal  flow  whereas  radial  stresses  increase 
significantly  due  to  combustion. 

Profiles  of  measured  axial  and  radial  fluctuations  along  the  stag¬ 
nation  plane  showed  a  good  degree  of  symmetry  with  respect  to 
the  stagnation  point  as  shown  in  Fig.  18.  Fluctuation  levels  in¬ 
creased  in  the  radial  direction  in  the  presence  of  combustion.  By 
contrast,  isothermal  flows  [6]  only  showed  an  increase  in  the  axial 
Reynolds  stress  component  from  0.20 14  at  the  stagnation  point  to 
0.29L4  at  a  radial  distance  of  approximately  0.5D,  while  the  radial 
component  decreased  slightly  from  0.12 Ub  to  0.10L4.  The  increase 
of  the  radial  stress  component  for  combusting  flows  is  consistent 
with  the  increased  radial  movement  illustrated  by  the  broader 
probability  density  function  for  the  stagnation  point  location 
shown  in  Fig.  11. 

Goh  et  al.  [22]  outlined  the  extension  of  the  POD  technique  to 
combusting  flows  by  applying  the  density  segregation  technique 
to  provide  conditional  POD  data.  Results  for  the  current  cases  are 


shown  in  Figs.  19  and  20.  As  shown  in  Fig.  19,  the  introduction 
of  impact  plates  in  combination  with  conventional  turbulence  gen¬ 
erating  plates  leads  to  a  substantial  decrease  in  turbulence  intensi¬ 
ties.  Hence,  such  a  combination  is  not  feasible  when  elevated 
turbulence  levels  are  sought.  By  contrast,  the  introduction  of  im¬ 
pact  plates  has  only  a  modest  impact  on  the  generated  turbulence 
levels  when  fractal  grids  are  used.  The  implication  is  that  the  influ¬ 
ence  of  the  flow  conditions  upstream  of  the  fractal  grid  has  a  much 
reduced  impact.  It  is  further  readily  apparent  that  the  uncondi¬ 
tional  POD  results  are  biased  towards  the  reactants  and  that  the 
presence  of  the  flame  reduces  the  levels  of  turbulence  in  the  reac¬ 
tants  as  compared  to  the  isothermal  flow.  Results  obtained  for  six 
methane  and  propane  flames  are  shown  in  Fig.  20  for  fractal  grids 
combined  with  impact  plates  and  the  trends  are  consistent  with 
the  methane  case  discussed  above. 

3.2.3.  Reaction  progress  variable  statistics 

The  first  two  moments  of  the  reaction  progress  variable  statis¬ 
tics  derived  from  PIV-images  obtained  using  aluminium  oxide 
seeding  in  combination  with  the  flame  front  detection  algorithm 
are  shown  in  Fig.  21.  The  methane  and  propane  flames  featured 
Ub  =  4.0  m/s  and  ethylene  flames  Ub  =  8.0  m/s  with  equivalence  ra¬ 
tios  between  0.70  and  0.90.  The  average  flame  front  location  is  evi¬ 
dent  from  the  mean  reaction  progress  variable.  It  can  be  seen  that 
the  flames  move  closer  together  when  approaching  the  lean 
extinction  limits.  Leaner  fuel/air  mixtures,  e.g.,  methane  or 
propane  at  1 f>  =  0.80  or  ethylene  at  < f>  =  0.70  show  a  reduction  of 
the  mean  progress  variable  at  the  nominal  stagnation  point. 
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Fig.  14.  Comparison  of  velocity  fluctuations  for  isothermal  (a)  and  reacting  (o) 
flows  along  the  burner  centreline  for  methane/air  flames  with  </>  =  0.9  and  a  bulk 
velocity  of  4.0  m/s.  Top:  Perforated  plates  without  impact  plates.  Middle:  Perforated 
plates  with  impact  plates.  Bottom:  Fractal  grids  with  impact  plates. 


Fig.  16.  Mean  normalised  axial  velocities  (left)  along  the  burner  centreline  and 
mean  normalised  radial  velocities  (right)  at  the  stagnation  plane  for  methane  (top) 
and  propane  (middle)  flames  at  a  bulk  velocity  of  4.0  m/s.  Ethylene  flames  (bottom) 
are  shown  for  a  bulk  velocity  of  8.0  m/s.  Equivalence  ratios  of  0  =  0.9,  o  =  0.8  and 
A  =  0.7. 


Centreline  [x/D]  Centreline  [x/D] 

Fig.  15.  Effect  of  an  increasing  a  circular  filtering  threshold  on  centreline  profiles  of 
the  normalised  axial  mean  velocity  (upper),  axial  (middle)  and  radial  (lower) 
velocity  fluctuations  for  a  methane/air  flame  with  </>  =  0.8  (right)  and  </>  =  0.9  (left)  at 
a  bulk  velocity  of  4.0  m/s.  Solid  line:  No  threshold.  Dashed-dotted  line  filtering  with 
a  3.5  mm  threshold.  Dashed  line  filtering  with  a  2.5  mm  threshold.  All  profiles  were 
obtained  using  fractal  grids  with  fitted  impact  plates. 


Fig.  17.  Comparison  of  the  normalised  axial  (left)  and  radial  (right)  velocity 
fluctuations  along  the  burner  axis  for  methane  (o)  and  propane  (a)  flames  at  a  bulk 
velocity  of  4.0  m/s.  Top  row:  </>  =  0.9.  Middle  row  </>  =  0.8.  Bottom  row:  Normalised 
axial  (left)  and  radial  (right)  velocity  fluctuations  along  the  burner  axis  for  ethylene 
flames  at  </>  =  0.8  (0)  and  </>  =  0.7  (□)  at  8.0  m/s. 
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Fig.  18.  Normalised  axial  (left)  and  radial  (right)  velocity  fluctuations  for  methane 
(top)  and  propane  (middle)  flames  along  the  stagnation  plane  for  a  4.0  m/s  bulk 
velocity.  Results  for  ethylene  flames  (bottom)  were  obtained  for  a  bulk  velocity  of 
8.0  m/s.  Stoichiometries  of  0  =  0.9,  o  =  0.8  and  a  =  0.7. 


O) 

a3 


as 

rs 

E 

IS 

O 


Mode  number  [ 


Mode  number  [ 


Fig.  19.  Cumulative  energy  in  modes  via  conditional  and  unconditional  POD, 
normalised  and  multiplied  by  the  mean  turbulence  energy  in  the  flow  field  for 
methane-air  flames  with  0  =  0.90,  averaged  over  all  images  in  each  set.  The 
turbulence  energy  was  for  comparison  purposes  defined  to  be  (1  )(u'u'  +  i/i/). 
Conventional  turbulence  generating  plates  with  (TGP-1P)  and  without  (TGP-NOIP) 
fitted  impact  plate.  Fractal  grid  with  (FRAC-IP)  and  without  (FRAC-NOIP)  fitted 
impact  plate.  Dotted  line:  Isothermal  flow;  Dashed-dotted  line:  Reactants;  Dashed 
Line:  Products;  Solid  line:  Unconditional  POD. 


Accordingly,  reactants  can  also  be  found  along  the  stagnation 
plane.  In  general,  propane  flames  appear  more  wrinkled  and 
curved  towards  their  reactant  streams  while  methane  flames  are 
flatter,  thinner  and  with  the  flame  brushes  closer  together  (cf.  Luff 
et  al.  [17]  and  Yoshida  et  al.  [65]).  Variations  in  the  separation 
distance,  also  found  by  Sato  [66]  and  Ishizuka  and  Law  [67]  for 


Fig.  20.  Cumulative  energy  in  modes  via  conditional  and  unconditional  POD, 
normalised  and  multiplied  by  the  mean  turbulence  energy  in  the  flow  field 
averaged  over  all  images  in  each  set.  The  turbulence  energy  was  for  comparison 
purposes  defined  to  be  (l)(u'u'  +  i/v1).  Fractal  grids  with  a  fitted  impact  plate  used 
in  all  cases.  Left  column:  Methane-air  flame  with  0  =  0.80  (top),  0.90  (middle)  and 
1.0  (bottom).  Right  column:  Propane-air  flame  with  0  =  0.80  (top),  0.90  (middle) 
and  1.0  (bottom).  Dotted  line:  Isothermal  flow;  Dashed-dotted  line:  Reactants; 
Dashed  Line:  Products;  Solid  line:  Unconditional  POD. 


laminar  flames,  is  likely  due  to  differences  in  thermochemistry 
and  transport  properties  with  the  impact  of  the  Lewis  number  dis¬ 
cussed  below. 

A  comparison  of  methane  and  propane  flames  shows  that  for 
the  same  bulk  strain  rate,  the  twin  methane  flames  were  on  aver¬ 
age  closer  together  with  a  higher  probability  of  finding  reactants  at 
the  nominal  stagnation  point.  It  should,  however,  also  be  noted 
that  the  propane  flames  extinguished  at  an  equivalence  ratio  of 
0.79  whereas  methane  flames  were  stable  even  below  0.75.  At 
equal  stoichiometries,  propane  has  a  higher  laminar  burning  veloc¬ 
ity  and  a  higher  adiabatic  flame  temperature.  The  higher  Lewis 
number  (see  Table  1)  may  be  the  reason  for  the  earlier  extinction 
of  the  propane  flames.  Ethylene  flames  withstood  bulk  velocities 
in  excess  of  8.0  m/s  with  progress  variable  statistics  showing  sim¬ 
ilar  trends  to  the  other  fuels  when  subjected  to  a  much  larger  total 
strain  as  also  shown  in  Fig.  21.  The  overall  behaviour  of  the  fuels 
suggests  that  both  chemical  rate  and  transport  processes  affect 
progress  variable  statistics  at  turbulent  Reynolds  numbers  >300. 


3.2.4.  Conditional  velocities  and  scalar  fluxes 

The  reaction  progress  variable  information  obtained  from  the 
flame  front  detection  algorithm  was  used  to  generate  conditional 
statistics  following  a  procedure  similar  to  that  by  Kalt  et  al.  [68], 
Hence,  the  mixture  was  treated  as  reactants  if  the  instantaneous 
average  value  of  the  reaction  progress  variable  in  an  aligned  inter¬ 
rogation  window  was  less  than  50%  and  otherwise  as  products.  The 
associated  spatial  distribution  of  velocity  vectors  in  combination 
with  the  progress  variable  profiles  made  it  possible  to  obtain  con¬ 
ditional  velocities  and  scalar  fluxes.  The  corresponding  laminar 
flame  properties,  such  as  reactant  (pr)  and  product  (pp)  densities 
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Fig.  21.  Comparisons  of  the  mean  progress  variable  (c)  (left)  and  its  variance  (cT7) 
(right)  for  methane  (top)  and  propane  (middle)  flames  with  a  4.0  m/s  bulk  velocity. 
Ethylene  flames  (bottom)  at  a  bulk  velocity  of  8.0  m/s.  Stoichiometries  of  </>  =  0.9 
(0),  </>  =  0.8  (o)  and  </>  =  0.7  (a). 


and  Lewis  numbers  for  the  deficient  reactant  were  computed  using 
a  subset  of  a  detailed  chemical  kinetic  mechanism  [45]  and  are 
presented  in  Table  1. 

Statistics  were  calculated  along  the  stagnation  point  streamline 
for  the  conditional  axial  velocities  of  reactants  (ur)  and  products 
(Op),  as  well  as  the  scalar  flux  (i Fc7).  The  data  was  analysed  in  the 
context  of  the  Bray-Moss-Libby  formalism  [69]  with  the  mean  slip 
velocity  (us).  defined  as  the  difference  between  the  mean  product 
and  reactant  velocities  [70],  also  calculated.  Results  for  ur,up  and 
us  obtained  for  </>  =  0.80  and  0.90  are  presented  in  Fig.  22  for  meth¬ 
ane  flames  and  in  Fig.  23  for  propane  flames.  The  quantities  are 
plotted  against  the  corresponding  mean  reaction  progress  variable. 
The  slip  velocity  is  almost  constant  in  the  region  0.2  ^  c  ^  0.7  for 
all  cases.  Towards  the  leading  edge  of  the  flames  the  reactants 
decelerate  leading  to  a  rapid  decrease  in  the  slip  velocities.  A  decel¬ 
eration  of  both  reactant  and  product  velocities  is  apparent  further 
towards  the  product  side.  The  results  suggest  that  the  flow  along 
the  centreline  can  be  classified  in  three  regions.  Two  of  the  regions 
(c  ^  0.2  and  c  ^  0.7)  show  deceleration  and  an  increase  in  pres¬ 
sure  while  in  a  narrow  region  in  physical  space,  corresponding  to 


Fig.  22.  Conditional  reactant  and  product  velocities  (left)  and  slip  velocities  (right) 
for  methane  flames,  (o)  Reactant  (ur)  and  (□)  product  (up)  velocities  for  the  upper 
nozzle,  (x)  Reactant  (ur)  and  (+)  product  (up)  velocities  for  the  lower  nozzle.  Mean 
slip  velocities  (us)  for  the  upper  (a)  and  lower  (*)  nozzles.  Top  row:  </>  =  0.8.  Bottom 
row:  </>  =  0.9. 


Fig.  23.  Conditional  reactant  and  product  velocities  (left)  and  slip  velocities  (right) 
for  propane  flames,  (o)  Reactant  ( ur )  and  (□)  product  (up)  velocities  for  the  upper 
nozzle,  (x)  Reactant  (ur)  and  (+)  product  (up)  velocities  for  the  lower  nozzle.  Mean 
slip  velocities  (us)  for  the  upper  (a)  and  lower  (*)  nozzles.  Top  row:  </>  =  0.8.  Bottom 
row:  </>  =  0.9. 


Table  1 

Values  of  density  used  for  calculations  of  the  Favre  averaged  terms.  Reactant  ( pr )  and  product  ( pp )  densities.  The  Lewis  number  was  defined  for  the  deficient  reactant  (fuel)  as 
0LeD  =  Gscfr'jGpr  and  laminar  burning  velocities  obtained  from  calculations  using  detailed  chemistry  [45].  Rates  of  dissipation  in  the  reactant  streams  (er)  were  obtained  using  the 
conditional  dissipation  method  defined  by  Goh  et  al.  [22]  as  derived  from  the  work  by  George  and  Hussein  [30]. 


Fuel 

</> 

pr  (kg  m  3) 

Pp  (kg  m  3) 

vr  x  10  5  (m2  s  ’) 

er  (m2  s-3) 

Si(m2  s->) 

si*,,  (-) 

ch4 

0.80 

1.134 

0.177 

1.60 

49.9 

0.262 

0.96 

ch4 

0.90 

1.129 

0.168 

1.59 

53.0 

0.325 

0.96 

ch4 

1.0 

1.124 

0.162 

1.59 

44.9 

0.368 

0.96 

c3h8 

0.80 

1.192 

0.176 

1.50 

48.9 

0.293 

1.82 

c3h8 

0.90 

1.193 

0.168 

1.49 

45.7 

0.362 

1.82 

c3h8 

1.0 

1.194 

0.164 

1.48 

42.5 

0.412 

1.80 

c2h4 

0.70 

1.178 

0.179 

1.53 

197 

0.309 

1.32 

c2h4 

0.80 

1.178 

0.164 

1.52 

210 

0.420 

1.31 

2432 


K.H.H.  Coh  et  al.  / Combustion  and  Flame  161  (2014)  2419-2434 


roughly  3  mm  (0.2  ^  c  <  0.7),  the  dominant  effect  is  a  flow  accel¬ 
eration  caused  by  the  heat  release  across  the  flame  brush. 

The  turbulent  scalar  fluxes  along  the  stagnation  point  stream¬ 
line  are  presented  as  a  function  of  the  reaction  progress  variable 
in  Fig.  24  and  axial  distance  in  Fig.  25.  The  latter  figure  also  shows 
the  estimated  total  uncertainties  and  it  is  apparent  that  the  scatter 
increases  significantly  for  the  stoichiometric  flames  and  hence  the 
corresponding  conditional  velocity  analysis  is  subject  to  increased 
uncertainties.  However,  it  is  apparent  that  the  flames  are  increas¬ 
ingly  in  the  transition  region  between  gradient  (GD)  and 
counter-gradient  (CGD)  diffusion  as  the  equivalence  ratio  is 
increased  and  the  experimental  scatter  may  be  influenced  by  the 
delicate  balance  of  the  turbulent  transport  mechanisms.  The 
increased  scatter  is  also  partly  due  to  more  prevalent  preheating 
(^50  K)  of  the  reactant  streams,  particularly  from  the  upper  nozzle, 
as  a  result  of  the  increased  rate  of  heat  release.  The  impact  of  the 
latter  is  also  evident  for  the  ethylene  flames.  However,  the  change 
in  the  sign  of  the  scalar  flux  can  readily  be  analysed  using  BML 
theory  [69]  and  is  accordingly  determined  by  the  slip  velocity  as 
shown  in  Eq.  (6). 

u"c"  =  c(l  -  c)(up  -  ur)  =  c(l  -  c)us  (6) 

A  negative  scalar  flux  suggests  gradient  diffusion  and  a  positive 
value  counter  gradient  diffusion  and  the  current  experimental  res¬ 
olution  is  sufficient  to  show  the  narrow  region  of  GD  at  the  leading 
edge  (cf.  Fig.  24).  Veynante  et  al.  [71  ]  derived  the  Bray  number  as  a 
criterion  for  the  transition  between  GD  and  CGD  in  freely  propa¬ 
gating  flames  with  subsequent  refinements  provided  [72,73],  Chen 
and  Bilger  [74]  further  analysed  the  transition  from  GD  to  CGD 
using  Bunsen  flames.  Lindstedt  and  Vaos  [3]  derived  an  alternative 
relationship  based  on  the  balance  between  mean  strain  (produc¬ 
tion  term  tensor  P,j)  and  mean  pressure  gradient  effects  (tensor 
&ij).  However,  in  the  turbulent  opposed  jet  configuration  the 
geometry  imposed  pressure  gradient  exerts  a  significant  influence 
on  the  transition  from  GD  to  CGD.  Efforts  have  been  made  to  ac¬ 
count  for  such  effects  [68,72],  though  a  comprehensive  evaluation 
of  the  resulting  model  formulations  present  difficulties  in  the  cur¬ 
rent  context.  However,  the  current  data,  see  Fig.  25,  shows  a  grad¬ 
ual  transition  between  GD  and  CGD  regimes  as  cf>  is  increased.  The 
result  is  consistent  with  increases  in  the  heat  release  gradually 
overcoming  the  adverse  effect  of  the  pressure  gradient  leading  to 
an  expanding  local  region  of  CGD  with  increasingly  positive  fluxes. 


0.0  0.5  1.0  0.0  0.5  1.0 

C  [-1  c  [-1 

Fig.  24.  Axial  scalar  fluxes  Wcr  for  methane  (left)  and  propane  (right)  flames  shown 
in  reaction  progress  variable  (c)  space,  o  Upper  nozzle.  □  Lower  nozzle.  Top  row: 
0  =  0.8.  Bottom  row:  ij>  =  0.9. 


Fig.  25.  Axial  scalar  fluxes  ire7  for  methane  (left)  and  propane  (right)  flames  shown 
in  physical  space.  Stoichiometries  of  0  =  0.8  (top),  0  =  0.9  (middle)  and  0  =  1.0 
(bottom). 


3.2.5.  Turbulent  burning  velocity  determinations 

Multiple  definitions  of  the  turbulent  burning  velocity  (e.g.,  Bray 
[70])  have  been  advanced  and  Driscoll  [75]  showed  that  definitions 
(e.g.,  [76,77])  are  not  equivalent.  The  matter  is  explored  quantita¬ 
tively  below.  Bray  [70  j  derived  the  turbulent  burning  velocity  (sT ,i ) 
as  twice  the  magnitude  of  the  slip  velocity  at  the  leading  edge  of 
the  flame.  Due  to  the  steep  gradients  in  slip  velocities  close  to 
the  leading  edge  the  direct  evaluation  can  result  in  significant 
uncertainties  while  the  determination  of  the  relevant  root  to  Eq. 
(7)  using  linear  interpolation  remains  amenable. 

Sr.i  =  -2  (us)M 

.-.  solve  us :  prur  +  2pus  =  0  (7) 

Lawn  and  Schefer  [77]  provided  a  simple  definition  of  the  tur¬ 
bulent  burning  velocity  (s72)  based  on  the  mean  flow  velocity  at 
the  leading  edge  (c  =  0.02)  as  shown  in  Eq.  (8).  From  a  theoretical 
point  of  view,  Catlin  and  Lindstedt  [78]  have  shown  that  the 
burning  velocity  eigenvalue  (zt)  is  strongly  dependent  upon  the  c 
iso-contour  (c*)  used  for  the  determination  with  A  =  0.54  for 
c*  =  10^,21  =  1.25  for  c*  =  l(T2,zt  =  1.63  for  c*  =  1CT3  and 
A  =  1.80  for  c*  =  1CT4  and  the  method  can  be  expected  to  have  a 
tendency  to  produce  lower  values. 

Sr, 2  =  M  Ic- =0.02  (8) 

Driscoll  [75]  proposed  the  use  of  the  local  minimum  in  the 
mean  velocity  profile  as  an  indication  of  the  leading  edge  location. 
However,  the  axial  mean  velocity  profiles  shown  in  Fig.  16  suggest 
that  local  minima  are  only  prominent  at  higher  equivalence  ratios. 
It  is  accordingly  proposed  here  that  the  local  minimum  in  the 
gradient  of  the  axial  velocity  profile  should  used  to  indicate  the 
location  of  the  leading  edge  of  the  flame  and  hence  used  to 
determine  the  burning  velocity  (s73)  as  shown  in  Eq.  (9).  The  veloc¬ 
ity  gradients  were  determined  by  backward  differencing  with  qua¬ 
dratic  least  square  fits  used  for  noise  reduction. 
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Table  2 

Values  of  different  turbulent  burning  velocities  and  mean  axial  RMS  fluctuations  (u') 
evaluated  as  averages  of  upper  and  lower  nozzles.  Values  of  sTt\ ,  sTi2  and  Sr,3  are 
derived  from  Eqs.  (7)-(9).  Uncertainties  in  u'  were  evaluated  as  unbiased  standard 
deviation  from  the  mean  of  corresponding  localised  values  for  sT,i,sT<2  and  Sr,3.  The 
theoretical  values  for  sT/u’  =  ( sL  +st,kpp)/u'  were  derived  using  Eq.  (10)  with  values 
for  sL  and  er  obtained  from  Table  1  with  the  reactant  value  of  vK  =  (vrer)1/4.  The  value 
of  CR  was  obtained  as  4.0 /<rLeD  [83]  with  resulting  values  for  sT/u'  given  below. 


Fuel 

0 

u!  (m/s) 

sT/u'  (- 

)  sT,,/u'  (-) 

St, i/u'  (-) 

ST.3/U'  (-) 

ch4 

0.80 

0.90  ±  0.03 

2.45 

2.32  ±0.17 

2.36  ±0.16 

2.72  ±  0.09 

ch4 

0.90 

0.92  ±  0.01 

2.75 

2.87  ±0.16 

2.45  ±  0.07 

2.78  ±  0.05 

ch4 

1.00 

0.90  ±  0.02 

3.01 

3.17  ±0.07 

2.63  ±  0.20 

2.90  ±  0.33 

C3H8 

0.80 

0.89  ±  0.02 

2.01 

2.17  ±0.17 

2.31  ±0.01 

2.61  ±0.15 

C3H8 

0.90 

0.89  ±  0.02 

2.29 

2.79  ±  0.04 

2.52  ±  0.08 

2.79  ±  0.04 

c3h8 

1.00 

0.89  ±  0.02 

2.50 

2.78  ±  0.24 

2.74  ±0.10 

2.96  ±0.14 

c2h4 

0.70 

1.49  ±0.00 

1.90 

2.31  ±0.01 

2.40  ±  0.03 

2.75  ±  0.07 

c2h4 

0.80 

1.49  ±0.04 

2.25 

2.48  ±0.12 

2.38  ±  0.02 

2.98  ±  0.06 

St, 3  =  (9) 

dx 

The  centreline  profiles  of  p,c,u,ur,us  and  dii/dx  were  used  to 
determine  turbulent  burning  velocities  through  the  alternative 
definitions  (st,i,St,2  and  S73)  and  statistics  from  the  lower  nozzle. 
The  results  are  presented  in  Table  2  and  values  agree  to  within 
20%  -  comparatively  close  given  the  relatively  wide  spread  high¬ 
lighted  in  the  reviews  by  Driscoll  [75]  and  Lipatnikov  and  Chomiak 
[79  .  Overall,  the  values  of  Sj-3  are  larger  than  the  corresponding  sT2 
determination  (based  on  c*  =  0.02)  as  the  point  of  evaluation  was 
consistently  closer  to  the  true  leading  edge  with  mean  velocities 
corresponding  larger.  However,  it  can  be  seen  that  there  are  excep¬ 
tions  with  values  of  Sj-.i  lower  than  sT,2  at  an  equivalence  ratio  of 
0.80,  which  may  reflect  the  evaluation  difficulties  associated  with 
the  leading  edge  slip  velocity  as  discussed  above.  For  comparison 
purposes  [22],  a  theoretical  estimate  can  readily  be  provided  via 
a  KPP  analysis  [80,81]  applied  to  a  fractal  flame  surface  area  based 
rate  expression  [3,82],  The  resulting  scaling  ( AKPP  =  2)  is  given  in 
Eq.  (10)  with  the  required  parameters  listed  in  Table  1. 

St, kpp  =  Akpp \J Cr(?> /2)(CI1/ Vk )  if'  VK  =  (vrer)'^4  (10) 

In  the  theoretical  estimation  of  the  turbulent  burning  velocity 
(st  =  sL  +sTKPp),  the  velocity  fluctuations  (tT)  normal  to  the  flame 
propagation  direction  were  used.  The  values  for  sT/u'  shown  in 
Table  2  provide  reasonable  agreement  with  a  tendency  to  under¬ 
predict  the  experimental  values.  The  latter  is  particularly  notice¬ 
able  for  ethylene.  Fractal  grids  have  been  found  to  substantially 
enhance  turbulent  diffusivities  [26]  due  to  the  multi-scale  nature 
of  the  generated  turbulence  and  a  modest  reduction  in  a‘Sc  would 
improve  the  overall  agreement. 

4.  Conclusions 

The  approach  to  extinction  of  fuel  lean  turbulent  premixed 
flames  in  fractal  grid  generated  multi-scale  turbulence  was  studied 
for  methane,  propane  and  ethylene  flames  using  P1V  with  the  flows 
to  both  nozzles  seeded  with  1  pm  silicon  oil  droplets  or  3  pm  A1203 
particles.  The  obtained  data  include  velocity  and  reaction  progress 
variable  statistics  as  well  as  conditional  velocities  and  scalar  fluxes. 
The  fractal  grids  were  used  to  ameliorate  limitations  associated 
with  conventional  perforated  plate  generators  with  the  turbulent 
Reynolds  number  range  moved  from  50-120  to  130-318.  Further¬ 
more,  optimal  grid  configurations  were  determined  with  particular 
emphasis  on  reducing  the  impact  of  the  flow  upstream  of  the 
turbulence  generators  in  order  to  facilitate  simpler  boundary  con¬ 
ditions  for  computational  studies.  It  has  been  shown  that  the  cur¬ 
rent  fractal  grids  have  the  notable  advantage  of  removing  the 


sensitivity  to  upstream  conditions  while  maintaining  significantly 
increased  levels  of  turbulence. 

Dilatation  due  to  the  heat  release  across  the  flame  front  induces 
a  positive  velocity  gradient  which  leads  to  an  acceleration  of  the 
flow.  The  reduction  in  pressure  gradient  due  to  dilatation  results 
in  a  drop  in  the  axial  fluctuations  along  the  axis  for  all  measured 
fuel  types.  Leaner  flames  showed  less  acceleration  across  the  flame 
front  due  to  lower  heat  release.  The  same  trends  were  found  for  the 
mean  radial  velocity  and  the  fluctuations  along  the  stagnation 
plane.  Radial  strain  rates  along  the  stagnation  plane  reached  their 
maximum  value  at  a  radial  distance  of  0.3D  for  all  fuel  types.  The 
radial  strain  rates  measured  in  the  reacting  flow  show  an  increase 
~100%,  compared  to  isothermal  measurements  [6],  that  is  directly 
related  to  the  higher  radial  velocity  gradients  caused  by  heat 
release  along  the  stagnation  plane. 

The  resulting  flow  structures  were  further  analysed  using  POD 
and  CPOD  with  results  confirming  that  conventional  (uncondi¬ 
tional)  POD  is  biased  towards  the  reactants.  The  use  of  CPOD  en¬ 
ables  a  more  detailed  analysis  of  the  flow  field  with  the  impact 
of  the  flame  on  both  reactant  and  product  fields  readily  segregated. 
The  technique  shows  that  conventional  POD  leads  to  a  hybrid 
biased  towards  the  reactants.  Probability  density  functions  were 
determined  for  the  instantaneous  location  of  the  stagnation  point 
to  eliminate  the  possibility  of  low  frequency  bulk  motion  distorting 
velocity  statistics.  The  PDFs  of  flame  curvature  were  determined 
using  a  multi-step  flame  front  detection  algorithm  with  estimates 
of  the  turbulent  burning  velocity  provided  along  with  a  discussion 
of  alternative  determination  methods.  The  data  sets  show  that 
fractal  grids  generate  multi-scale  broadband  turbulence  and  pres¬ 
ent  an  opportunity  for  a  systematic  evaluation  of  calculation  meth¬ 
ods  for  premixed  turbulent  flames  that  undergo  a  transition  from 
non-gradient  to  gradient  turbulent  transport  while  approaching 
extinction.  It  is  further  shown  that  reasonable  agreement  with 
models  for  the  turbulent  burning  velocity  is  obtained  for  the  cur¬ 
rent  set  of  flames.  The  lean  extinction  limits  for  methane  and  pro¬ 
pane  flames  suggest  an  influence  of  the  Lewis  number  and  were, 
respectively,  at  </>  =  0.73  and  c/>  =  0.79  for  a  fixed  bulk  velocity  of 
4.0  m/s  with  ethylene  flames  extinguishing  at  </)  =  0.66  at  8.0  m/s. 
It  is  expected  that  the  current  data  sets  will  be  of  particular  value 
in  the  development  of  generalised  models  that  account  for  the 
transition  between  gradient  and  non-gradient  transport  in  the 
presence  of  imposed  pressure  gradients. 
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